Theory of magnetization plateaux in the Shastry-Sutherland model 
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Using perturbative continuous unitary transformations, we determine the long-range interactions 
between triplets in the Shastry-Sutherland model, and we show that an unexpected structure de- 
velops at low magnetization with plateaux progressively appearing at 2/9, 1/6, 1/9 and 2/15 upon 
increasing the inter-dimer coupling. A critical comparison with previous approaches is included. 
Implications for the compound SrCu2(BC>3)2 are also discussed: we reproduce the magnetization 
profile around localized triplets revealed by NMR, we predict the presence of a 1/6 plateau, and 
we suggest that residual interactions beyond the Shastry-Sutherland model are responsible for the 
other plateaux below 1/3. 

PACS numbers: 05.30.Jp, 03.75.Kk, 03.75.Lm, 03.75.Hh 



O 

(N 



i 

CO 

03 



i 

C 

O 
o 



> 
o 
m 

o 

OO 

o 



Following the discovery of magnetization plateaux in 
the layered copper oxide SrCu2(B03)2[l, 2], a lot of ac- 
tivity has been devoted to the properties in a magnetic 
field of the 2D spin- 1/2 Heisenberg model known as the 
Shastry-Sutherland modcl[3] and defined by the Hamil- 
tonian: 



H 



J' S t -S J + J Si-Sj-B^S* 



<h0> 



where the <C i,j 3> bonds build an array of orthogonal 
dimers while the < i,j > bonds are best seen as inter- 
dimer couplings (see Fig. 1). For J'/J smaller than a 
critical ratio of order 0.7, the ground state of the model 
is exactly given by the product of dimer singlets, and the 
magnetization process can be described in terms of po- 
larized triplets It 1 ) = | ||) on the dimers interacting and 
moving on an effective square lattice [4, 5]. These triplets 
can be described as hard-core bosons, and the magneti- 
zation plateaux correspond to Mott insulating phases. 

All theoretical approaches agree on the presence of 
magnetization plateaux at 1/3 and l/2[4, 5, 6, 7, 8], in 
agreement with experiments[l, 9]. However, the struc- 
ture below 1/3 is rather controversial. On the exper- 
imental side, the original pulsed field data have only 
detected two anomalies interpreted as plateaux at 1/8 
and 1/4[1], but the presence of additional phase tran- 
sitions and of a broken translational symmetry above 
the 1/8 plateau has been established by recent torque 
and NMR measurements up to 31 T[10, 11]. The pos- 
sibility of additional plateaux has been pointed out by 
Sebastian et al[9], who have interpreted their high-field 
torque measurements as evidence for plateaux at 1 j q with 
2 < q < 9 and at 2/9. On the theoretical side, the sit- 
uation is not settled either. The finite clusters available 
to exact diagonalizations prevent reliable predictions for 
high-commensurability plateaux, and the accuracy of the 
Chcrn-Simons mean-field approach initiated by Misguich 
et al.[7] and recently used by Sebastian et al.[9] to ex- 
plain additional plateaux is hard to assess. The essen- 
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FIG. 1: Shastry-Sutherland lattice and definition of the 2- 
body interactions. V n is the coefficient of the 2-body interac- 
tions between the thick dimer and the dimer labeled V„ . 



tial difficulty lies in the fact that, since plateaux come 
from repulsive interactions between triplets, an accurate 
determination of the low-density, high-commensurability 
plateaux requires a precise knowledge of the long-range 
part of the interaction, which could not be determined 
so far. 

In this Letter, we combine perturbative continuous 
unitary transformations (PCUTs)[12] with an analysis 
of the effective hard-core boson model reformulated as a 
spin model to investigate the magnetization process of 
the Shastry-Sutherland model. As we shall see, below 
1/3, the results are well controlled up to quite large val- 
ues of J' j J, providing compelling evidence in favour of a 
rich and unexpected plateau structure at small magneti- 
zation. 

The PCUT transforms the Shastry-Sutherland model 
into an effective model conserving the number of elemen- 
tary triplets (triplons[13]). The relevant processes for 
the physics in a finite magnetic field have maximum to- 
tal spin and total S z . Other spin channels relevant for 
spectroscopic observables have been studied earlier [14]. 
The general form of the effective Hamiltonian obtained 
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by the PCUT reads therefore: 

n— 2,4,6--- ri,"- ,r n 

where the r/s are sites of the square lattice formed by the 
J-bonds, while the hard-core boson operator b\ creates a 
polarized triplet l^ 1 } at site r. The coefficients CVi,r 3 ,— )Tn 
are obtained as a series in J'/ J. We have kept all terms 
with up to 3 creation and annihilation operators and all 
4-body interactions (n ri n r2 n r3 n ri ) that first appear at 
order < 8. For the 2-body interactions (n ri n r2 ), which 
dominate the physics at low densities (see below), we 
keep more terms, namely those that first appear at order 
less or equal to 10. The coefficients (7 riir2i ... >Tn are evalu- 
ated up to order 15 for the 2-body interactions and up to 
order 12 for the other terms, and they are then extrapo- 
lated using Pade or DlogPade extrapolants. The result- 
ing Hamiltonian contains more than 15000 processes. 

Clearly it is not possible to describe all the terms of 
-f/eff. To give an idea of its structure, the dominant 
terms of each type are sketched in Fig. 2. In the range 
< J'/ J < 0.5, the relative importance of the various 
terms (as deduced from their coefficients or from their 
contribution to the ground state energy as estimated be- 
low) follows simple trends. The pure density-density in- 
teraction terms have by far the largest coefficients. So 
one expects the magnetization curve to consist mainly 
of plateaux separated by jumps, or by very narrow in- 
termediate phases. Secondly, the magnitude of all terms 
decreases when the distance between sites involved in the 
interaction increases. As a consequence, at low density, 
interaction terms involving more than two particles give a 
negligible contribution to the energy, which is completely 
dominated by 2-body interactions. Thirdly, the standard 
two-site hopping appears only at order 6[15] and is really 
small, a consequence of the frustration of the model, and 
the kinetic part of the Hamiltonian is dominated by the 
correlated hopping (n ri 6j 9 6 r3 ), which allows a particle to 
hop only if there is another particle nearby [4, 14]. Such 
a correlated hopping has been recently shown to strongly 
favour supersolid phases [16], and supersolid rather than 
supcrfluid intermediate phases are possible. 

Given their importance at low density, we show in 
Fig. 3 the evolution with J' / J of the 2-body interactions 
defined in Fig. 1. Although at low J' /J the interactions 
smaller than V4 may be neglected, this is not true for 
larger J'/ J where the higher order terms V3, V5, and 
Vj (appearing at order 6) become important and con- 
tribute to the formation of low-density plateaux. For 
these terms, the bare series and the Pade extrapolations 
are basically indistinguishable below J'/ J = 0.5. Beyond 
that value and up to J' / J ~ 0.63, various Pade extrap- 
olations still give consistent results for these 2-body in- 
teractions. Remarkably, we find that local interactions 
involving more than two particles become very strong 
for J' /J > 0.6 close to the phase transition and are very 
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FIG. 2: Largest process of each type in the Hamiltonian H c g 
and its amplitude G ru ... , r „ in units of J at J' / J = 0.5. 
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FIG. 3: Coefficients of the extrapolated 2-body interactions 
as a function of J /J. Inset: Different extrapolants (solid 
lines) as well as the bare series (dashed lines) for V3 and Vs. 

hard to extrapolate. We therefore restrict the discussion 
to J' I J < 0.5 where the expansion is well controlled. 

The effective Hamiltonian H c r has both positive and 
negative off-diagonal terms, which prevents its investiga- 
tion by Quantum Monte Carlo (QMC). So, to determine 
its phase diagram, we map it onto a spin 1/2 model us- 
ing the Matsubara-Matsuda representation[17] of hard- 
core bosons S + = b. S~ = b', S z = 1/2 — b'b, and we 
use a classical approximation where the spins are treated 
as classical vectors of length 1/2. Given the complexity 
of the model, the energy is then minimized numerically 
on finite size clusters with periodic boundary conditions. 
To allow for high commensurabilities of various symme- 
tries, all non-equivalent clusters with up to 32 sites for all 
J'/ J (64 for J'/ J = 0.5) have been tested and compared. 
The possible phases are identified using symmetry argu- 
ments. A state is supcrfluid if it breaks the U(l) gauge 
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FIG. 4: Well converged plateaux as a function of fi and J'/ J. 
The boson density n is equal to the magnetization in units of 
the saturation value, and the chemical potential fi is equal to 
the magnetic field B. 



symmetry [z axis spin rotation), it is solid if it breaks the 
translationnal symmetry (S z is non-uniform) , and it is a 
supersolid if it breaks both symmetries. At the classical 
level, the relevant order parameters are the in-plane com- 
ponent of the spins (y '(S x ) 2 + (Sv) 2 = \{b)\) for the su- 
perfluid and the static structure factor for the solid. This 
classical approximation has been shown to be remark- 
ably good in cases where its predictions could be com- 
pared with QMC[16, 18]. More precisely, this approach 
slightly overestimates the tendency to form plateaux but 
plateaux not present in this method are very unlikely to 
occur. 

The resulting phase diagram is shown in Fig. 4, and 
the structure of the plateaux in Fig. 5. The solid lines 
in the phase diagram denote results that are fully con- 
verged with respect to the terms kept in the Hamiltonian 
(see Fig. 6). Well converged domains are then connected 
by dashed lines. The phase diagram is dominated by a 
series of plateaux, at 1/3 and 1/2 (not shown) already at 
very small J'/ J, then by plateaux at 2/9, 1/6, 1/9 and 
2/15. Whether these plateaux are separated by jumps or 
intermediate phases (with possibly additional plateaux) 
cannot be decided on the basis of the present calculation. 

The actual magnetization profile inside the various 
plateaux is depicted in Fig. 5. The spin density profiles 
are obtained by transforming the relevant observable in 
the PCUT formalism [12]. In all cases, the building brick 
is a triplet and its two up-down neighboring dimers, in 
agreement with the interpretation of Cu NMR in the first 
plateau of SrCu 2 (B0 3 ) 2 [2]. 

Let us now compare the present results with previ- 
ous works. Momoi and Totsuka used pertubation theory 
to third order and logically found only plateaux at 1/3 
and 1/2 [4]. Miyahara and Ueda used a phenomenologi- 
cal form of the long-range 2-body interactions in a model 
without kinetic energy to successfully determine possible 
structures inside various plateaux with emphasis on the 
plateaux at 1/8, 1/4 and 1/3 reported in SrCu2 (1303)2 [6]. 
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FIG. 5: Spin density (S z ) profile of the main plateaux at 
J'/ J = 0.5. Full (empty) circles corresponds to magnetization 
along (opposite to) the magnetic field. The radius of the 
circles is proportional to the magnetization amplitude. The 
blue line shows the unit cell compatible with the periodicity 
of the state. For the 2/15 plateau, two structures have the 
same energy within the error bars of the method. 
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FIG. 6: Magnetization curve at J'/ J — 0.5: comparison of 
the results obtained by keeping terms with up to 2 creators 
(dotted blue curve), 3 creators (solid blue curve) and 4 cre- 
ators (black curve). Well converged plateaux used in Fig. 4 
are indicated explicitly. Inset: Energy of the 1/9, rhomboid 
1/8, 2/15, and 1/6 plateaux as a function of /i at J'/J = 0.5. 
The error bars from the DlogPade extrapolation are smaller 
than the line width. 
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However, their model exhibited several other plateaux, 
and the approach was not set up to be predictive re- 
garding the actual plateaux stabilized in the model. Fi- 
nally, a Chern-Simons theory where spins are mapped 
onto fermions has been developed by Misguich et al.[7], 
who only found plateaux at 1/4, 1/3 and 1/2, and by Se- 
bastian et al. [9] , who allowed for non- uniform mean-field 
solutions and found plateaux at 1/q with q = 2, 9 and 
at 2/9. There are obvious similarities between the results 
of Sebastian et al and ours. In particular, all our plateaux 
except the 2/15 are present in their series, and this excep- 
tion could be due to the fact that this plateau requires 
a unit cell not considered in Ref. 9. But some aspects 
of their results, for instance a well developed plateau at 
1/5, are definitely ruled out by our analysis and must be 
considered as artefacts of the method. So further work 
is needed to understand the reliability and pitfalls of the 
Chern-Simons approach to magnetization plateaux. 

Finally, let us discuss the implications for 
SrCu2 (B03)2. The most surprising aspect of our 
results is the absence of plateaux at 1/8 and 1/4, which 
have been observed in SrCu2(B03)2- This is especially 
so for the 1/8 plateau, which has not only been observed 
with pulsed field, but has been extensively studied in 
steady fields with torque[ll] and NMR[2]. To convince 
the reader that this plateau is indeed absent in the 
Shastry-Sutherland model at J' / J = 0.5, we have 
plotted in the inset of Fig. 6 the energy of various 
low-density plateaux including the 1/8 plateau with 
rhomboid unit cell [8]. Clearly, the 1/8 structure is never 
stabilized. 

Taking for granted that the experimental identification 
of these 1/8 and 1/4 plateaux is correct (a point which 
might actually deserve further investigation in view of 
the rather different absolute scales reported in pulsed and 
steady magnetic field experiments), we can think of two 
possible origins of this discrepancy. The first one is that 
the physics changes dramatically when J' / J approaches 
the critical value where the dimer product wave-function 
is no longer the ground state but is replaced by a phase 
with plaqucttc order[19]. This seems unlikely however in 
view of the first order nature of this transition. 

The other possibility is that the plateau structure is in- 
fluenced by residual couplings beyond the simple SU(2) 
Shastry-Sutherland model. There arc a priori three types 
of residual couplings, all of the same order of magni- 
tude (J/100): in plane Dyzaloshinskii-Moriya interac- 
tions (both inter- and intra-dimer), further neighbor in 
plane exchange couplings, and inter-plane exchange cou- 
plings. In view of the very large parameter space gener- 
ated by these couplings, a systematic study of their in- 
fluence will be difficult. Recent ab-initio results might be 
helpful in restricting the relevant portion of this parame- 
ter space though [20]. This is left for future investigation. 
Here it is interesting to notice that the 1/4 and rhomboid 
1/8 plateau[8] share with the 1/6 plateau the peculiarity 



to have all triplets parallel to each other, and to maxi- 
mize the number of pairs of up-down dimers which are 
nearest neighbors on the same sublattice. Hence, if resid- 
ual interactions indeed favor the 1/8 (resp. 1/4) plateaux 
over the (1/9,2/15) pair (resp. 2/9), they are expected 
to further stabilize the 1/6 plateau, which thus appears 
as a very strong candidate for an additional plateau in 
SrCu 2 (B0 3 )2. 

To summarize, a systematic investigation of the effec- 
tive Hamiltonian describing the triplets induced by an 
external field in the Shastry-Sutherland model has lead 
to the prediction of an unexpected series of additional 
plateaux below 1/3 at 2/9, 1/6, 1/9 and 2/15 that ap- 
pear progressively upon increasing the inter-dimcr cou- 
pling. High-order perturbation theory was essential in 
identifying the correct sequence of plateaux, a remark 
that opens interesting perspectives for other models. Re- 
garding SrCu2(B03)2, our results suggest to have a closer 
look at the magnetization below 1/3, with two issues in 
mind: the absolute value of the magnetization at the 
plateaux currently assumed to be at 1/4 and 1/8, and 
the presence of a 1/6 plateau. 
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